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Rare particle decays that are either forbidden within 
the standard model of particle physics (SM), or are ex- 
pected to have very small branching ratios provide excel- 
lent signatures with which to look for new physics and 
allow to probe subatomic processes that are beyond the 
reach of direct searches. The decays B®^ — > e + y~ [l[ 
are forbidden within the SM, in which lepton number 
and lepton flavor are conserved. However the observa- 
tion of neutrino oscillations indicates that lepton flavor 
is not conserved. To date, no lepton flavor violating 
(LFV) decays in the charged sector such as B9 ■> — ► e + yT 
have been observed. These decays are allowed in mod- 
els where the SM has been extended by heavy singlet 
Dirac neutrinos The LFV decays are also allowed 
in some physics scenarios beyond the SM, such as the 
Pati-Salam model Q and supersymmetry (SUSY) mod- 
els (ij. The grand-unification theory by J. Pati and A. 
Salam predicts a new interaction to mediate transitions 
between leptons and quarks via exchange of spin-1 gauge 
bosons, which are called Pati-Salam leptoquarks (LQ), 
that carry both color and lepton quantum numbers [3|. 
The lepton and quark components of the leptoquarks are 
not necessarily from the same generation |5j, [6J , and the 
decays B® — > e + and B° — > e + pT can be mediated 
by different types of leptoquarks. Processes involving 
flavor-changing neutral currents (FCNCs) can occur in 
the SM only through higher-order Feynman diagrams 
where new physics contributions can provide a signifi- 
cant enhancement. Compared to B9s — > fi + [i~ 0], the 

FCNC decays of B9^ — > e + e~ are further suppressed by 
the square of the ratio of the electron and muon masses 
(rrie/m^) 2 . The SM expectations for branching ratios of 
B9 n — > e + e~ arc of the order of 10 
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In this Letter we report on a search for the LFV de- 
cays B9. — ► e+/i~ and the FCNC decays B^ — ► e + e _ , 

using a data sample corresponding to 2 fb _1 of integrated 
luminosity collected in pp collisions at y/s = 1.96 TeV. 
With no evidence for either the LFV or FCNC decays, we 
set upper limits on their branching ratios using the com- 
mon reference decay B° — > A' + 7r~, which has a precisely- 
known branching ratio. This is the first time a search for 
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B® — > e + e~ has been performed. 

A detailed description of the CDF II detector can be 
found in Ref. 0] . Here we give a brief description of the 
detector elements most relevant to this analysis. Charged 
particle tracking is provided by a silicon microstrip de- 
tector together with the surrounding open-cell wire drift 
chamber (COT), both immersed in a 1.4 T axial mag- 
netic field. The tracking system provides precise ver- 
tex and momentum measurement for charged particles 
in the pseudorapidity range < 1.0 [ToT ]. Surround- 
ing the tracking system are electromagnetic (CEM) and 
hadronic sampling calorimeters, arranged in a projective 
geometry Drift chambers and scintillation counters are 
located behind the calorimeters to detect muons within 
\r)\ < 0.6 (CMU) and 0.6 < \rj\ < 1.0 (CMX). 

We use a data sample enriched in two-body -B-decays 
selected by a three-level trigger system using the ex- 
tremely fast tracker [IH at level-1, and the silicon ver- 
tex trigger [l2| at level- 2. The trigger requires two 
oppositely-charged tracks, each with a transverse mo- 
mentum pt > 2 GcV/c, and an impact parameter (l3j 
0.1 < do < 1 mm. It also requires the scalar sum of 
the transverse momenta of the two tracks to be greater 
than 5.5 GeV/c, the difference in the azimuthal angles 
of the tracks 20° < Aip < 135°, and a transverse decay 
length [14j L xy > 200 /urn. At the level-3 trigger stage, 
and in the offline analysis, the trigger selections are en- 
forced with a more accurate determination of the same 
quantities. In the off-line analysis, additionally we re- 
quire: the B-meson isolation I > 0.675 [l5j|, the pointing 
angle A0 < 6.3° [l(|, and a tighter selection of L xy > 375 
/mi. These three thresholds were optimized in an unbi- 
ased way to obtain the best sensitivity for the searches 
using the procedure described in Ref. [17|. 

Electron and muon identification is applied in the selec- 
tion of B9 s j — > e + y~ and-B° s ) — > e + e~ decay modes. The 
electron identification [3] requires that both the specific 
ionization (dE / dx) measured in the COT, and the trans- 
verse and longitudinal shower shape as measured in the 
CEM, be consistent with the hypothesis that the particle 
is an electron. The performance of electron identification 
is optimized using pure electron samples reconstructed 
from 7 — > e + e~ conversions and hadron and muon sam- 
ples from D° — > A'~7r + , A — * pir~ , and J/tp — * 
decays. We find the identification efficiency to be around 
70% for electrons. The muon identification starts from 
tracks in the COT that are extrapolated into the muon 
detectors and are required to match hits in the muon 
systems. The muon selection is fully efficient for muons 
with pt > 2 GcV/c in CMU or CMX. 

The mass resolution a m of fully-reconstructed B— 
meson decays to two charged particles is about 
28 McV/c 2 . Energy loss due to bremsstrahlung by elec- 
trons generates a tail on the low side of the mass distri- 
bution. This tail is more prominent for the B9 % — > e + e _ 
channels, where two electrons are involved. We define 
search windows of (5.262-5.477) GeV/c 2 for B° s -► e+fx' 
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and (5.171-5.387) GeV/c 2 for B° -► e+/z~- These 
correspond to a window around the nominal values of 
the B® and B° masses [l9| of approximately ±3<r m . 
To recover some of the acceptance loss due to electron 
bremsstrahlung for the B9 s — > e + e~ channels, we choose 
wider and asymmetric search windows ranging from 6 
a m below to 3 a m above the nominal values of the 
Bg and B a masses. The search windows arc (5.154- 
5.477) GeV/c 2 for the B° s and (5.064-5.387) GeV/c 2 for 
the B°. The sideband regions (4.800-5.028) GeV/c 2 and 
(5.549-5.800) GeV/c 2 are used to estimate the combina- 
torial backgrounds. 

The background contributions considered include com- 
binations of random track pairs and partial B decays that 
accidentally meet the selection requirement (combinato- 
rial) , and hadronic two-body B decays in which both final 
particles are misidentified as leptons. The combinatorial 
background is evaluated by extrapolating the normalized 
number of events found in the sidebands to the signal 
region. The double-lepton misidentification rate is deter- 
mined by applying electron and muon misidentification 
probabilities to the number of two-body decays found in 
the search window. 

Figure [1] shows the invariant mass distribution for 
candidates. We observe one event in the B® mass 
window, and two events in the B mass window, con- 
sistent with the estimated total background of 0.8 ± 0.6 
events in the B® search window, and 0.9 ± 0.6 in the B° 
window. The combinatorial background in both chan- 
nels is estimated to be 0.7 ± 0.6 events. The number of 
events where two tracks are misidentified as electron and 
muon is estimated to be 0.09 ± 0.02 for the B® case and 
0.22 ± 0.04 for the B° case. 

Figure [5] shows the invariant mass distributions for 
e + e _ candidate pairs where both tracks were identified 
as electrons. We observe one event in the B® mass win- 
dow, and two events in the B° mass window. We esti- 
mate the total background contributions to be 2.7 ± 1.8 
events in both the B® and B° mass windows. The domi- 
nant contribution comes from combinatorial background: 
2.7 ± 1.8 compared to the contribution where both tracks 
are misidentified as electrons: 0.038 ± 0.008 for both B® 
or B°. 

We use the reference decay B° — > K + 7r~ to set a limit 
on B(B® — ■> e + £~) (where t is either e or /i), using the 
following expression: 



e + n 



B(B° S 



-£-\ - N ( B °- 



) 



■e+l-)-g(B Q ^if + Q-/ d // 3 



■N(B°-t K+tt-) 



The expression for the B° channels is identical, except 
that the ratio of fc-quark fragmentation probabilities: 
fd/fs is n °t present. In the expression, N(B° S — > e + £~) 
is the calculated upper limit on the number of B® — > 
e + i~ events, N(B° — ► K + ir~) is the observed num- 
ber of events from the reference channel B° — > K + ir~ , 
B(B° -> K+n-) = (19.4 ±0.6) x 10~ 6 is the branch- 



ing ratio for the B° — > K + tt decay, and e B e l +e _ is 
the detector acceptance and event selection efficiency for 
reconstructing B® — > e + l~ decays relative to that for 
B° -> K+tt-. The value of f d /f s is 3.86 ± 0.59, where 
the (anti-)correlation between the uncertainties has been 
accounted for [2(| ■ To calculate the detector acceptance, 
we use simulated events with a detailed simulation of 
the CDF II detector and event selection. We obtain 



el 



_ = 0.207 ±0.016, egb'^ e+ = 0.210 ± 0.012, 



0.129±0.011, and e^iLe+e- = 0.128±0.011. 



B°- 
rel 

e B°^e+e~ 

The uncertainties listed above are the combined statisti- 
cal and systematic uncertainties. The later include un- 
certainties from detector fiducial coverage, electron and 
muon identification efficiencies, detector material deter- 
mination, B? s s px spectrum, and £?P s lifetimes. The 
reference channel B° — > K + ir~ has been reconstructed 
using the same selection criteria except lepton identifica- 
tion. We find 6387 ± 214 B° -> K+tt- events, using a 



21] 



fitting procedure similar to that described in Ref. 

The upper limit on the branching ratio in each search 
window is obtained using the Bayesian approach [3] , as- 
suming a flat prior, and incorporating Gaussian uncer- 
tainties into the limit. The total systematic uncertain- 
ties, listed in Table [H are used as input for the limit cal- 
culation. Table |TT] lists the upper limits we obtain on the 
branching ratios at 90% (95%) credibility level (C.L.). 

Within the Pati-Salam leptoquark model, the follow- 



ing relationship between the B(B? 



leptoquark mass (Mlq) can be derived [B[ 



[i ) and the 



B(B° M - eV) = ™ 2 (Af LQ ) 



F 2 B 3 m 3 B , 



T B a 



where R 



m 

00 



y a s (rn t ) J \a. 



(™t) 



. The val- 

m b y a 3 (m t ) J \a s (m b ) ) 

ues and uncertainties of the quantities used in the calcula- 
tion of Mlq are the following (ljj : the top-quark mass mt 
(171.2 ± 2.1 GeV/c 2 ), the bottom quark mass m b (4.20 
± 0.17 GeV/c 2 ), the charm quark mass m c (1.27 ± 0.11 
GeV/c 2 ), the B°-meson mass m B o (5.27953 ± 0.00033 
GeV/c 2 ), the £°-meson mass m B o (5.3663 ± 0.0006 
GeV/c 2 ), the B°-meson lifetime t b o (1.530 ± 0.009 ps), 
the 5°-mcson lifetime r B a (1.470 ± 0.027 ps), the cou- 
pling strength Fro (0.178 ± 0.014 GcV), and F B o (0.200 
± 0.014 GeV)|22j. For the strong coupling constant we 
use a s (Mz°) = 0.115, which is evolved to Mlq using the 
Marciano approximation [23| assuming no colored parti- 
cles exist with masses between mt and Mlq- Using the 
limits on the branching ratios listed in Table HB we calcu- 
late limits on the masses of the corresponding Pati-Salam 
leptoquarks of M L q(B s -► e+pr) > 47.8 (44.9) TeV/c 2 
and M L q(B° e + n~) > 59.3 (56.3) TeV/c 2 at 90 
(95)% C.L. Figure [3] shows the limit and the relation 
between the leptoquark mass and the branching ratio for 
the B2 meson. 
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TABLE I: Values used to calculate the limits on B(B? S \ — > e+ /i ) and B(B? S \ — > e+e ) and their uncertainties. 



Source 


Values 


13 (Bg -> e+/i- 


) B{B U -» e+M" 


) -> e+e- 


-) B{B» -> e+e-) 


N(B {> -► /V+7T-) 


6387 ±214 


3.4% 


3.4% 


3.4% 


3.4% 




(19.4 ±0.6) x 10" 


6 3.1% 


3.1% 


3.1% 


3.1% 


/b°//b° 


3.86 ±0.59 


15.3% 




15.3% 






0.207 ±0.016 
0.210± 0.012 


7.6% 


5.9% 








0.129 ±0.011 






8.9% 






0.128 ±0.011 








8.9% 


Total 




17.7% 


7.5% 


18.3% 


10.0% 



> 



<U 
> 
LU 



_ Sideband 



+ 

CD 

T 

o 
CD 



+ 
CD 

T 

o tn 
CO 



Sideband 



3.8 4.9 5.0 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 

M(e + |T) GeV/c 2 

FIG. 1: Invariant mass distribution of e+ [iT pairs for events 
where one track passed the electron identification and the 
other track the muon identification. The B® (B°) search win- 
dow is indicated by the solid (dotted) line. The sideband 
regions are indicated by the dashed lines. 







TABLE II: Branching ratio limits at 90(95) % C.L. 

BjB^ e + fi~) < 2.0 (2.6) x HF 7 
B{B° -> e+fj, ) < 6.4 (7.9) x 10~ 8 
B(B° S -> e+e-) < 2.8 x 10" 7 
B(B° -> e+e") < 8.3 x lO" 8 



B(B^->en) <2.0(2.6)x10~ 7 90 (95)% C.L. 
M LQ (B° -> e n) > 47.8(44.9) TeV/c 2 90 (95)% C.L. 

G. Valencia, S. Willenbrock, PRD 50:6843 (1994) 
with: T(BJ = 1 .470 ± 0.027 ps 

m(B s ) = 5.3663 ± 0.0006 GeVIc* 
n\ ' =4.20 ±0.17GeV/c 2 
= 0.200 ± 0.014 GeV 




M LQ (B"->en) [TeV/c^] 



FIG. 3: Leptoquark mass limit corresponding to the 90 (95) 
% C.L. on B(B° a -> e + n~). The error band is obtained by 
varying the values entering the theoretical calculation within 
their uncertainties. The uncertainties stemming from approx- 
imating a s are not included. 



Sideband 



CD 
+ 
CD 

t 
o 
CO 



CD 
+ 
CD 

T 

o tn 
CO 



Sideband 



5.1 5.2 5.3 5.4 5.5 

M(e + e) GeV/c 2 



FIG. 2: Invariant mass distributions of e+e~~ pairs for events 
where both tracks passed the electron identification. The B® 
{B°) search window is indicated by the solid (dotted) line. 
The sideband regions are indicated by dashed lines. 



In summary, we report on a search for the lepton flavor 
violating decays B^ — ► e+ and the flavor changing 
neutral current decays B, v — > e+e- using data corre- 
sponding to 2 fb^ 1 of integrated luminosity collected in 
pp collisions at y/s = 1.96 TeV. This is the first search 
for B° 



z+e~ decays. We observe no evidence for these 
decays and set limits that are the most stringent to date. 
These results represent a significant improvement com- 
pared to the previous measureme nt 12411 b y CDF and the 
best results from i?-Factories [25|, |26|. |27| . 
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